THEORY OF ADSORBATE INDUCED 
SURFACE RECONSTRUCTION ON 
W(IOO) 



Kari Kankaala *'^, Tapio Ala - Nissila'^'^'**, and See- Chen Ying-'- 

^ Physics Department, Box 1843, Brown University, 
Providence, RI 02912, USA 



o\ 

Oh! 
(D 

00 

^ ■ *Center for Scientific Computing, P.O. Box 40, 

SF-02101 Espoo, Finland 

' Department of Electrical Engineering, Tampere University 

O 

a^ 
o 

^ ■ ** Research Institute for Theoretical Physics, 



of Technology, P.O. Box 527, SF-33101 Tampere, Finland 



University of Helsinki, 



^ • Siltavuorenpenger 20 C, SF-00170 Helsinki, Finland 

a : 

O . Abstract 
O 

- , , , We report results of a theoretical study on an adsorbate induced surface 

X 

^ , reconstruction. Hydrogen adsorption on a VF(IOO) surface causes a switching 

transition in the symmetry of the displacements of the W atoms within the ordered 
c(2 X 2) phase. This transition is modeled by an effective Hamiltonian, where 
the hydrogen degrees of freedom are integrated out. Based on extensive Monte 
Carlo renormalisation group calculations we show that the switching transition 
is of second order at high temperatures and of first order at low temperatures. 
This behavior is qualitatively explained in terms of an XY model where there is 
an interplay between four and eight fold anisotropy fields. We also compare the 
calculated phase diagrams with a simple mean field theory. 
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1. Introduction 

Structural phase transitions have been observed on a number of surfaces. One of 
the most extensively studied systems undergoing a surface reconstruction is the 
clean VF(IOO) surface (see e.g. [1,2,3] and refences therein). It has a very rich 
phase diagram with a reconstructive phase transition from a (1 x 1) phase into a 
c(2x2) phase [1,4,5,6]. The microscopic structure of the low temperature phase was 
unraveled early by Debe and King [5] and it has since then been accepted that this 
phase involves alternate displacements of the W atoms along the < 11 > direction 
in the surface plane. A number of first principles total energy calculations [e.g. 
7,8,9] have verified the inherent instability of the ideal (1 x 1) phase towards this 
reconstruction, and the Debe-King structure is found to have the lowest energy. 
There was some dispute on the nature of the high temperature phase, and thus also 
on the details of the structural phase transition between these two phases. Earlier, 
there was both experimental [1,10,11,12] and theoretical [13] evidence which was 
interpreted in favor of an ordered high temperature phase where the tungsten 
atoms would vibrate about their ideal bulk positions. The other possibility 
was that the tungsten atoms would be displaced from their bulk positions also 
in the high temperature phase but the displacements would point into random 
directions, thus forming a disordered phase [7,14]. Recent experimental [6,15,16] 
and theoretical [17,18,19,20] studies have corroborated the disordered nature of 
the high temperature phase, and it is now generally accepted that the transition 
is of order-disorder type instead of a displacive one. 

While the clean surface of VF(IOO) has been studied extensively both experimen- 
tally and theoretically, the infiuence of adsorbates on the reconstruction has been 
studied to a far less extent and in most cases a thorough understanding of these 
phenomena is missing. Among the adsorbates, the adsorption of hydrogen has been 
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studied in most detail [1,4,10,22,23,26,27,44]. It was found that small amounts 
of hydrogen enhance the inherent instability of the clean surface, and increase 
the critical temperature of the order-disorder transition [21]. The structure with 
adsorbed hydrogen was initially believed to be the same as the clean surface re- 
constructed c(2 X 2) phase or very similar to it [1,5,22,23]. It was later shown 
[10] and verified [4,44] that the hydrogen adsorbed surface has a different struc- 
ture when the coverage exceeds 9 ~ 0.1. Namely, the clean c(2 x 2) has a p2mg 
symmetry whereas the hydrogen rich (coverage ^ > 0.1) c(2 x 2) — H phase has 
a c2mm symmetry. The symmetry change was cleverly verified by Griffiths et 
al. [24] in a low energy electron diffraction (LEED) spot intensity study. They 
studied the difference in the relative intensities of two half-order beams, namely 
(|, i) and (|, |) at a normal incidence. The difference between these two beams 
is finite in the reconstructed clean surface p2mg phase, and it vanishes when the 
entire surface has switched into the c2mm phase. This symmetry switching was 
initiated at a coverage of 9 ~ 0.04 and completed at 9 ~ 0.16. Similar symmetry 
switching has been observed in an infrared spectroscopy study [25]. The c2mm 
phase is interpreted as due to a structure in which the W atoms have alternate 
displacements along the < 10> or <01 > directions. 

The hydrogen induced switching of the VF(IOO) surface between these two 
structures is the main subject of the present work. Our study is motivated by 
theoretical [26] and experimental [24] indications that the switching transition 
could be of first order at low temperatures, instead of being an XY like continuous 
transition [27] . We will first briefiy discuss the theoretical model for this work. For 
the clean surface, we adopt a lattice dynamical Hamiltonian which has recently 
been shown to describe the critical properties of the VF(IOO) surface very accurately 
[17,18,19,20]. It corresponds to an XY model with an intrinsic cubic anisotropy 
field. To this model then, we add the short range interactions of hydrogen with its 
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nearest neighbor W atoms. The hydrogen degrees of freedom are then exphcitly 
integrated out leaving an effective Hamiltonian. This new Hamiltonian is shown 
to contain anisotropy fields of all orders allowed by symmetry, and we argue that 
the nature of the switching transition at low temperatures is dictated by the 
interplay between effective fourth and eighth order anisotropy fields. We then 
proceed to calculate the details of the phase diagram, using first simple mean 
field arguments to locate the line of switching transitions on the (T, /x) - plane. 
To study the nature of the transition in detail, we present results of extensive 
Monte Carlo Renormalisation Group calculations. These results show that at 
high temperatures the transition is continuous, while below a multicritical point 
Tt it becomes discontinuous. Finally, we summarize our results and discuss the 
relevance of our work to experimental studies. Preliminary results have been 
previously reported [28] . 

2. Theoretical Models 

As discussed above, it has been experimentally verified that hydrogen induces 
a switching transition where the H atoms adsorbed on bridge sites pin the 
displaced W atoms towards themselves resulting in displacement orientations 
along either < 10 > or < 01 > direction, instead of the clean surface < 11 > 
direction. This reconstruction has been subject to theoretical investigations {e.g. 
[3,8,9,26,27,29,30,31]) but to a far less extent than the clean W{100) surface. 

The models discussed below, as well as our model belong to the class of so called 
lattice dynamical models which are written in terms of lattice displacements. 
When studying the driving forces for the transition, microscopic details have to be 
included in the model. However, lattice dynamical models are usually sufficient 
to study the critical properties of the systems under consideration. The role of 
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the detailed electronic degrees of freedom in these models has been discussed by- 
various authors [3,32,33]. 

2.1 Lattice Dynamical Models 

To study the effect of an adsorbate on the switching transition on the W^(IOO) 
surface, Lau and Ying [27] suggested a lattice dynamical model Hamiltonian 

H = H°'^^° + H^*^. (1) 

In their work, the clean surface term H'^^®^" is of the type suggested by Fasolino 
et al. [34] , with the details left unspecified. The crucial point is that for the clean 
VF(IOO) surface, it has the symmetry of an XY model with cubic anisotropy. The 
intrinsic anisotropy of the clean surface is such that the W atoms are displaced 
along the < 11 > direction in the reconstructed c(2 x 2) phase. The adsorbate 
part H^'^ in Eq. (1) was chosen to be 

^'^=\Y^ Jvo'^i'^r + ni'v{R^ + Ui- Ri') (2) 

which describes the effects of hydrogen on the structural transition. In Eq. (2), n^/ 
is an occupation number, and Ji>j' is an interaction term between hydrogen atoms. 
The second summation describes the interaction between hydrogen and tungsten 
atoms as a sum of pair potentials v where Ri' is the coordinate of a hydrogen 
adsorption site and ^^'s denote the ideal bulk positions of the surface tungsten 
atoms. The addition of ttj describes the spontaneous displacements of the surface 
atoms. 

Lau and Ying [27] started from Eq. (2) and found that the adsorbate induces a 
cubic anisotropy field that opposes the intrinsic clean surface one. This adsorbate 
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induced anisotropy will then increase in magnitude with coverage leading to the 
change of sign of the total anisotropy field. Hence the switching occurs from 
the phase in which W atoms are displaced along the < 11 > directions to a 
phase in which they are displaced along the < 10 > directions, favored by the 
opposing cubic term. This conclusion was supported by electronic band structure 
calculations [29] where it was suggested that even at low hydrogen coverages it 
would be energetically favorable for the surface atoms to align along the < 10 > 
or < 01 > directions in the low temparature c(2 x 2) phase with hydrogen atoms 
on bridge sites. Fasolino et al. [35] suggested the switching to occur at a coverage 
of ^ ~ 0.044 i^)"^ where a is the lattice constant and ^ is the surface coherence 
length. This value for switching agrees with LEED studies [24] . 

Ying and Roelofs [3,36] developed the model further by choosing for the clean 
surface Hamiltonian 



where Ui is the displacement of the surface W atom at the i lattice site, and 
and Uiy are the X and y components of the displacement, respectively. A and 
B are coefficients of a simple double well potential, Ci is the nearest neighbor 
interaction term, and h% is the clean surface cubic anisotropy field. For the clean 
VF(IOO) surface h\ < ^ which aligns the displaced W atoms along the < 11 > 
directions. Summation i in Eq. (3) goes over all lattice sites. 

Using (2) and (3), Ying and Roelofs obtained an effective H®*^ to describe the 
hydrogen effects by integrating out the hydrogen degrees of freedom in the mean 




(3) 
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field approximation, and working with a constant coverage description. However, 
this description breaks down if the transition happens to be of first order as we 
shaU discuss below. Also, the approximation of retaining only the lowest order 
anisotropy fields is not justified for strong H — W coupling as addressed below. 
In a closely related computer study, Sugibayashi et al. [26] used an extended 
anisotropic XY model to describe the clean surface, with a two - body lattice 
gas model for the adsorbate part. They presented results from their computer 
simulations for the switching transition, but were not able to conclude the nature 
of the transition at low temperatures. Their results seem to indicate, however, 
that there is a phase coexistence region at low temperatures between the < 11 > 
and < 10 > phases, which would be strong evidence for a first order transition. 

2.2 New Model for H/W(100) 

To be able to study a possible first order regime, and to give a more accurate 
description of the low temperature phase we have undertaken simulations of a 
new lattice dynamical model. The clean surface contribution in our model is that 
described by Eq. (3) to which we have added an additional eight fold anisotropy 
field: 



TTclean _ / ^^,2 , .,4 , oi,0.,2 „,2 

i 

-\- C±(^UixUjx ~\~ U^x^kx ~l~ ^iy^jy ~^ ^iy^ky) 

+ /i°cos(80O}, (4) 

where 0j = axcta,n{uiy / Uix) is the displacement angle. To calculate the total 
surface - adsorbate partition function, we will make the simplifying assumption of 
neglecting direct H — H interactions. Thus, we can write 
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with H^*^ = Ylii' iT'i' i^i' — A*); where the summation i goes over aU lattice sites and i' 
over all bridge sites. As usual P = 1/(A;bT), n^/ is the actual hydrogen occupation 
number at a bridge site i', Vi' is the interaction potential between hydrogen and 
surface tungsten atoms, is the chemical potential, and Ui is the displacement 
vector of a tungsten atom. Next, we sum over all hydrogen degrees of freedom to 
obtain 

Z = TV{,^}exp[-/3(H^*T)], (6) 

where H^^ is an effective Hamiltonian, which can be obtained in a straightforward 
fashion as 



jjeff ^ jj^lean _ ^ ^^^^^ ^ exp[-/?(V-, - fi)]}. (7) 

i' 

The summation i' goes over all bridge sites. They can also be labelled as i' = {i, v) 
in which there are two bridge sites u associated with each lattice site i. The 
interaction Vi' is obtained by expanding a pairwise interaction potential to lowest 
order in the W atom displacements. This gives us Vi^^ = a^j ■ Uj where the 
components of a are given by am — (— q;,0), aii2 = (ct, 0), 0^21 = (0, —a), and 
CKi22 = (0,0;), and Uj = {ujx,Ujy)'^ , and a is an effective interaction parameter 
between hydrogen and tungsten atoms. The relative positions of H and W atoms 
as well as displacement vectors are illustrated in Fig. 1. 

The Hamiltonian (7) is written in terms of the chemical potential whereas the 
experiments are carried out at a constant coverage. The description of the 
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adsorbate in terms of the chemical potential rather than the coverage allows us 
to better investigate the possibility of both a first and a second order switching 
transitions. The coverage can be calculated and is given by the expression 



in which the summation z, v goes over all bridge sites and the average < > denotes 
an average over configurations. Since there are two bridge sites for each W , the 
maximum allowed coverage is normalized to the value of two. 

2.3 Effect of the Anisotropy Fields in the Model for H/W(100) 

We can predict qualitatively the behavior described by H*'^ if we express it 
entirely in terms of the displacement angles 0j for fixed amplitudes instead of the 
displacement vectors ttj. If we retain only the isotropic nearest neighbor coupling 
terms and the leading anisotropy fields in the on - site potential we then obtain a 
simplified Hamiltonian 



Hfy = K J] cos(0, - - /i4(T, /x) COS Act), + hsiT, ij) ^ cos 9>4>i. (9) 



The anisotropy fields h^iT,fi) = hl{T,fi) + hi and hsiT.^i) = h^iT, fi) + 
where superscripts and a refer to clean surface and adsorbate induced terms, 
respectively. This Hamiltonian is an XY model with both four and eight fold 
anisotropics. 

The form of B.^^ together with the dependence of and hs on T and fi has 
important consequences for the expected behavior of the H/W{100) system since 
renormalisation group analysis [37] tells us that the four fold field is relevant at 





9 



all temperatures T < Tc- For small values of fj,, h4{T, fj,) < which favors the 
orientation of the W atoms along the < 11 > directions. However, since h%{T,ii) 
is a monotonically increasing function of jj, (for a fixed T), there will be a point 
(//sw)?sw) where h4{fisw,Tsyf) = 0. Beyond this, h/^iT^jj) > 0, and the minimum 
energy state is now given by W atoms displaced along the < 10 > directions. 
This simple analysis predicts that the switching transition in the (T, /i) - plane 
is determined by the condition /i4(//sw) 2^sw) = 0, and it is an XY transition of 
Kosterlitz - Thouless type [37,38,39]. Physically, for the clean surface, the intrinsic 
four fold field hi < indicating a preference of the displacements for the W atoms 
along the < 11 > directions. As H is adsorbed, it sits on the bridge sites and the 
H — W interaction energy is minimized if locally the displacements of the nearest 
neighbor W atoms switch and point along the < 10 > or < 01 > directions. This 
local distortion has a finite range. Globally, the average /i4(T, //) is still negative 
for low hydrogen coverages. However, the effective /i4(T, fx) now increases with the 
coverage or the chemical potential, refiecting these local distortions around each 
H atom. 

To understand the nature of the switching transition in more detail, we must also 
consider the effect of the eight fold anisotropy field hs (T, /i) which is induced by 
the H — W interaction term in (7). It is well known that the eight fold field is 
relevant at low temperatures but above some Tt < Tc it becomes irrelevant [37]. 
When the eight fold field is negative, it favors orientations = n^, n = 0, ...,7, 
while for hs{T, fi) > the preferred orientations are cf) = f+n^, n = 0, ...,7. This 
is demonstrated in Fig. 2. Consider first a fixed, negative hs- The eight fold field 
favors the same directions which are also favored by both positive and negative 
h^iT, ii). Thus, when /i4(Tsw, /Usw) ^ from the negative side a relevant hs{T, /i) 
field keeps the system at the orientation (f) = j. However, as soon as the four fold 
field is positive and finite, it dictates again the orientation of the system which 
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suddenly switches from = ^ to = 0. Thus, the transition is expected to be of 
first order at low temperatures when the eight fold field is relevant and negative. 
However, the situation is changed above Tt where hg becomes irrelevant. Then, at 
^4(T'sw, A*sw) = the system is a pure XY model. The phase boundary thus has 
a multicritical point Tt where Tt and Tg^ coincide for some This is precisely 
the behavior for our choice of H*^^®^"^ as in (4) with a zero intrinsic eight fold field. 
The reason is that the hydrogen adsorbate induced four fold and eight fold fields 
are of the same sign, and thus compatible with each other. 

For the case of a positive h%{T,ji) which competes with the four fold field, more 
complicated behavior can result below Tt. We expect then the phase boundary 
separating the < 11 > and < 10 > phases to open up into two Ising - like 
transitions with an intermediate phase in between [40,41]. This is a subject of 
a separate study. 

It should also be noted here that this discussion of four and eight fold fields 
is strictly speaking only qualitative for i7/VF(100). As we discuss in the 
Appendix, higher order fields are comparable in magnitude with /i4 and hs for 
a realistic H — W coupling. Therefore, the following mean field and Monte Carlo 
Renormalisation Group analyses are based on the full effective Hamiltonian of Eq. 
(7). In addition, although the switching transition occurs at /i4(T, //) = 0, this four- 
fold field should be a renormalized effective field in which the short wavelength 
details have been averaged away. Thus, even if we are able to extract correctly the 
bare /i4 field in the Hamiltonian, it still would not allow us to locate the switching 
point exactly. This is one more reason why the numerical work is absolutely 
necessary. The discussion above, however, is valuable in that it nicely brings about 
the physical ideas behind the change in the order of the switching transition. 
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3. Mean Field Estimate and the Choice of Parameters for the Switching 
Transition 

Before setting out to an extensive numerical study of our effective Hamiltonian, it is 
useful to study the mean field solution for the switching transition. The mean field 
approximation is exact at T = and thus at least at low temperatures it should 
yield reasonable results. At higher temperatures, the temperature fluctuations 
affect the transition, making the switching occur at a lower coverage or chemical 
potential than the mean field theory predicts. 

The energy difference between the two phases that we denote hy (f) — ^ and (f) — 
{N is the number of tungsten atoms), is: 



AE = 



h''^(<^=t)-h^"(<^ = o) 



TT 



eff/ 



(10) 



where 



H^'^l^ = j) =2Nhl + Nhl - ksT log[l + exp(-/9(K,,(0 =j)- 



(11) 



and 



H^ff(0 = 0) =Nhl - /csT^^logfl + exp{-P{V,,i{4> = 0) - //))] 

+ N{^ + ^-2C,). (12) 



This can be simplified by noting that there are only four different kinds of energy 
changes related to the bridge sites (total of 2A'") that can take place: 
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AE =2Nh1 

N 



T^^^^°4 1 + J 

, ^,_r l + exp(-/^(-V2«-/^)) l 



1 ^ I- Xpi^ ; yll^ V t 


I')) 


1 + exp(/3/i) 




1 + exp(-/?(-V2a 




1 + exp(-/?(2Q! - 




1 + exp(-/?(\/2Q; - 




1 + exp(/5//) 




1 + exp(-/3(v^Q; - 




1 + exp(-/3(-2Q; - 





- y^sTiog : : :z) 7/::. (13) 



In order to study in more detail the low temperature behavior of the system, we 
set AE = and seek solutions to it. These solutions become exact at T = in 
the absence of fluctuations. 

At T = 0, we get AE/N = |//sw + 2/i4 + Q; = in the interval -2a < //gw < -V^a. 
The solution for ii^w is then /Xgw = 2(— 2/i4 — a) where //gw is the value for the 
chemical potential where the switching between < 11 > and < 10 > phases occurs. 
Thus, by fixing the values of a and hi we can compare the prediction of the mean 
field theory with the experimentally observed values of coverage for the transition. 
From the results we can also deduce a constraint between the four fold anisotropy 
field hi and the interaction parameter a as ^\h1\ < (2 — V2)a which restricts the 
choice of parameters, as will be discussed below. 

In our model, there is also a solution for a higher chemical potential which 
corresponds to the unphysical situation of switching back to the < 11 > phase 
at a higher coverage. This is due to our omission of direct interactions between 
hydrogen adatoms on adjacent bridge sites. A direct repulsion between the 
adatoms on these sites would eliminate this unphysical transition. Experimentally, 
many complicated structures [4,25] have been observed at high hydrogen coverages 
and these are outside the present scope of study. 
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We have worked with two different sets of parameters. The first set (set I) is 
very close to that of Yoshimori and coworkers [26] and was chosen so that direct 
comparasion could be made with their results. We are particularly interested in 
answering the question whether the transition is of first order at low temperatures. 
Their results indicates this possibility but the conclusions were unclear. A notable 
point in this set of parameters is that the clean surface four fold field is very 
weak compared to other similar studies [17,18,19,20] or total energy calculations 
[8,9]. 

The other set of parameters (set II) is more consistent with experimental work and 
total energy calculations. A Hamiltonian based on this set of parameters has been 
successfully used to explain many clean surface phenomena [17,18,19]. However, 
to prevent a switching into displacements along directions other than < 10 >, we 
need to set the eight fold field /ig of Eq. 4 to a negative value. 

In what follows, we work with dimensionless displacements u and temperature 
T defined as u = UgU, T = TgT where Tg and Ug are the scale factors. The 
numerical values for the first set of parameters are A = Au^/Tg = —2,B = 
Bui/Ts = 8.8, Ci = Cml/Ts = 1.5, = -0.1, a = 4.5, and = 0. The 
numerical values for set II are A — Aul/Tg = — 10,-B = Bu^/Tg = 40, and 
Ci = Ciul/Tg = 3.75, ^2 = -1-85, a = 17, and = -0.6. The displacement 
amplitude u = (4Ci — A)/(B + 8/2.4) — 1 for parameter set I and is set to the 
value 1 for the second set. The critical temperature for the clean surface {p = —00) 
is about Tc ~ 1.3 for the first parameter set and about — 2.3 for the other set. 
The experimentally observed values for Tc ^ 230K and u ~ 0.2A yield the scale 
factors Us — 0.2 and Tg ^ 175 for the first set and ~ 100 for the second set. 
For simplicity, from here on throughout the paper we will always use the scaled 
values for the parameters, e.g. T for T. 
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The choice of a = 4.5 in the first set of parameters follows of Roelofs and 
Ying [3]. In their approach, the effective Hamiltonian was obtained by keeping 
terms only to fourth order in the displacements of the W atoms. The various 
parameters including a were then determined by comparing theoretical results 
with the experimental data. As we show in the Appendix, the higher order terms 
are quantitatively important and cannot be neglected. Indeed, the estimate for 
the clean surface anisotropy in ref. [3] is about an order of magnitude smaller 
than the current first principles results. Thus, the value of a = 4.5 can also 
no longer be viewed a reliable estimate. In the second set of parameters, the 
anisotropy field is much stronger and close to the first principles result. As a 
consequence, the constraint (2 — ^/2)a > 4:\h1\ implies that a much larger value of 
a is required. Our choice of ck = 17 yields a switching coverage in agreement with 
the experimental observation. It is also consistent with the infrared vibrational 
spectroscopy data when the new larger value of the anisotropy field is taken into 
account. 

Finally, we should note that a smaller value of | 1 makes the transition more XY 
- like and thus more difficult to analyze numerically whereas a larger value would 
decrease the fiuctuations and result in more Ising like behavior in the vicinity of 
the transition with smaller finite size effects. Larger absolute values of the four 
fold field hi also hold the system preferrably in the < 11 > direction and thus 
would move the switching to occur at a higher chemical potential. However, 
the parameters within the model cannot be changed arbitrarily as discussed 
above, and the value of imposes restriction to the value of the interaction 
parameter a. Increasing means increasing a, and this leads the switching to 
occur approximately at a constant coverage over a large temperature range. The 
qualitative behavior of these two parameter sets is the same both in the (T, fx) and 
in {T,9) planes (latter not sketched here). 

15 



4. Numerical Results 

To study the nature of the switching transition for the fuU Hamiltonian of Eq. 
(7), we have carried out extensive Monte Carlo renormahzation group (MCRG) 
simulations [42] . By studying the cumulants of the order parameter, we have been 
able to determine the order of the transition both at low and high temperatures. 
These results are supported by our additional studies of the order parameter, its 
distribution function, coverage, and the observation of strong hysteresis at low 
temperatures, as we will eludicate below. 

In the simulations, we have used the standard Metropolis updating scheme. They 
were carried out using a nonconserved order parameter and Glauber dynamics. 
The amplitudes of the tungsten atom displacements Ui were held fixed as it 
has been shown recently [17,18] that displacement amplitude fluctuations arc not 
important for the critical behavior of this system. We studied mainly two system 
sizes, namely 64 x 64 and 32 x 32. Computational time involved with larger 
systems becomes rather formidable. The averages were computed over 50 000 
configurations for parameter set I, and over 20 000 for set II. 

In the MCRG studies, the lattice is divided into blocks of size L, and we then 
study the cumulants of the moments of the order parameter distribution function. 
For this, we used two order parameters, namely 

i 

<*^>=z^E(-i)^"^^^W, (15) 

i 

where u^i and the X and y components of the displacement vector at the 

i^^ lattice site, and ix and iy are the x and y coordinates of the i^^ lattice site in 
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a block of size L. The coefficient (— is a phase factor. < > denotes an 
average over configurations, and the summation goes over all lattice sites. In the 
< 11 > phase both (14) and (15) have finite values whereas in the transition to 
<10>(<01>) phase < Ux > {< Uy >) vanishes. There is no preference for one 
phase over the other. By studying the behavior of the cumulants of these average 
displacement components, we can probe the order of the switching transition. 

The block cumulants used are defined by 



where i refers to x and y. The variation of these two cumulants as a function 
of the block size L depicts a flow diagram analogous to that of a renormalisation 
group method. It can be shown [42] that these cumulants approach zero above 
Tc as the block size increases. Below Tc, both these cumulants tend to nonzero 
values Ul ^ U* = \,Vl ^V* — ^, and at a second order transition they tend 
to nontrivial values. 

We also monitored the combined order parameter of Eqs. (14) and (15): 



< ^ ^(-l)(^-+'''^«..«,i, (18) 

i 

For the ordered < 11 > phase this order parameter has a finite value, and in the 
< 10 > phase it vanishes. 

If the transition is of first order, one of the block cumulants is expected to 
flow towards the ordered state value of | at all chemical potentials. As the 
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transition occurs suddenly between two ordered phases, one of the order parameter 
components does not vanish at the transition { < Uy > for the < 01 > phase and 
< Ux > for the < 10 > phase) but abruptly changes from one ordered state value 
to the other. The other component vanishes suddenly at the transition { < Ux > 
for the < 01 > phase and < Uy > for the < 10 > phase). 

One should note that if the transition is very weakly of first order, it is very difficult 
to distinguish between it and a second order transition. This is true in our case in 
the vicinity of the multicritical point where the order of the transition changes. 

In the following, we will present numerical data only for one parameter set (set 
I, Figs. 3 and 4) but will show the phase digram for both parameter sets in Fig. 
5. The reason for this is that results for both sets are qualitatively very similar, 
and the quantitative differences can be seen in the respective phase diagrams. For 
parameter set I, the switching occurs at approximately 9 ~ 0.02 which is lower 
than the experimentally observed value. For set II, we reproduce the switching at 
approximately ^ ~ 0.1 which is in good agreement with experiments [24]. 

4.1. High Temperature MCRG Results 

Fig. 3 (a) shows the behavior for parameter set I at high temperatures (T = 1.4), 
and confirms our predictions for the second order transition. The switching 
transition in this case is from the < 11 > state into the < 01 > state, i.e. < Ux > 
vanishes. The cumulants Ulx oi Ux show behavior typical for a second order 
transition: in the < 11 > phase [jj, < —12) the cumulants approach the fixed point 
of Ulx ~ 0.67 and in the < 01 > phase (// > 11.5) the cumulants approach zero. 
We estimate the transition to happen in the region where chemical potential is 
— 12 < jjL < —11. Inside this range of chemical potential, the cumulants approach 
a nontrivial fixed point (// ~ —11.5). The problematic behavior of the cumulants 
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at = — 14 is attributed to large fluctuations at this temperature and the vicinity 
of two phase transitions (order - disorder transition and switching). 

When studying the cumulants of < tt^ >, we see that in the ordered phase they 
approach the value of ULy ~ 0.67, and then near the transition they flow to a 
nontrivial flxed point at ji e [—12,-11.5]. At larger values of ji the cumulants 
again approach the ordered state value. The errors have been taken from the 
results from the largest lattice sizes. 

Based on these results, we conclude that the switching transition at Tg^ = 1.4 
occurs at /Xgw € [—12, —11.5], and it is of second order. Data for the parameter set 
II are very similar, and result in a switching transition for T = 2.4 (< Tc, high 
temperature) at fisw £ [—28.4, —28.0]. 

4.2 Low Temperature MCRG Results 

In Fig. 3 (b) we depict the behavior of the cumulants at Tg^ = 0.5 again for 
the parameter set I. At chemical potential values /j, < —8.8, the cumulants of 
Ulx all approach the ordered phase value. At /j, > —8.8, the cumulants vanish. 
The cumulant behavior of U Ly shows no change in behavior when passing through 
the transition but the cumulants approach the ordered state value of 0.67 at all 
chemical potentials. We interpret this as evidence of a flrst order phase transition 
into the < 10 > phase at this temperature. The strange behavior of Uj^y at 
fj, = —8.8 is probably to be due to strong metastable effects combined with slow 
dynamics of the system at these low temperatures. 

Based on similar arguments, the switching for set II occurs for T = 0.6 at 
yUsw G [—26.7, —26.4]. In this case the metastability effects seem to be less severe, 
and the accuracy of simulations is thus better even at relatively low temperatures. 
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4.3 Results from Other Quantities 

Father qualitative support for the MCRG results can be found in our coverage 
studies. In Fig. 4 the coverage as in Eq. (8) is plotted as a function of temperature 
for different values of the chemical potential. If the transition were of first order a 
discontinuity in the real coverage would be observed at the transition. However, 
as we work with a finite system, we observe only a very sharp rise in coverage at 
the transition. 

In the case of a second order transition, we would expect the coverage to increase 
gradually when the transition region is passed. The data in Fig. 4 is in qualitative 
agreement with the prediction of a first order transition at low temperatures 
and a second order transition a higher temperatures. The abrupt increase in 
coverage is very pronounced at T = 0.2. Similar behavior was also observed for 
the order parameter. The jump in coverage was very clear in simulations with set 
II parameters in the first order regime, and more gradual at higher temperatures 
reconfirming the qualitative similarity in simulations for both parameter sets. 

In addition, we have observed strong hysteresis and phase coexistence (as 
determined from the order parameter distribution function) at low temperatures 
to further support the scenario for a first order transition. Hysteresis studies have 
not been performed for the second parameter set. 

4.4 The Phase Diagram 

Based on the mean field theory and our numerical results, we sketch phase 
diagrams for both parameter sets in Figs. 5 (a) and 5 (b) for the line of switching 
transitions as induced by hydrogen adsorption on the VF(IOO) surface. The 
temperatures for the (1x1) — > c(2 x 2) transition both for the clean surface 
{jj, = — oo) and for the low hydrogen coverages have been determined from MCRG 
studies. The solid order - disorder line is only a guide to the eye. 
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In both Figs. 5 (a) and 5 (b), the dashed hne denotes the mean field solution. 
The solid line is the part of the switching transition where it is expected to be of 
second order, and the dotted line depicts the first order switching transition. The 
agreement with the mean field line is very good at low temperatures where the role 
of fluctuations is not important. At higher temperatures the deviations are due 
to temperature fluctuations which cause the system to undergo a phase transition 
at lower chemical potentials than our mean fleld results would imply. The mean 
fleld result for the parameter set II at higher temperatures deviates particularly 
strongly from the MCRG data. 

The whereabouts of the multicritical point where the transition changes its 
order could not be accurately pinpointed by our simulations. The locations of 
the multicritical points depicted in Fig. 5 are only schematic. We expect it to 
be somewhere between 0.5 < Tt < 1.0 for set I, and 1.0 < Tt < 1.5 for set II. 
The error bars shown depict the variance of numerical results from a series of 
simulation runs. 

5. Summary and Discussion 

To summarize, we have developed a model Hamiltonian to describe the adsorbate 
induced effects on the Vr(lOO) surface. The clean surface part consists of a lattice 
dynamical Hamiltonian, which describes the (1 x 1) — > c(2 x 2) transition as 
studied by Han and Ying [17,18]. To include the hydrogen induced effects, we have 
considered a simple model of interactions between W and H atoms, and integrated 
out the hydrogen degrees of freedom to obtain an effective Hamiltonian. This 
leads to an XY model with anisotropy fields of all orders allowed by symmetry. 
In particular, the interplay between the fourth and eighth order anisotropy fields 
indicates that the switching transition should be first order at low temperatures 
and of second order at high temperatures. This prediction was confirmed by Monte 
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Carlo renormalisation group calculations, and further corroborated by studies of 
other quantities. We also used the MCRG studies together with a simple mean 
field theory to map out the line of switching transitions on the (T, jj,) plane. 

We have carried out simulations with two rather different sets of parameters, and 
obtained qualitatively similar behavior for both sets. With set I, we have been able 
to reproduce results similar to those of Sugibayashi et al, and re - interpret their 
findings in terms of first and second order transitions. The transition coverage 
is somewhat lower than experimantally observed. The second parameter set is 
more in agreement with total energy calculations, and the results are qualitatively 
same as for the first set. These both parameter sets reveal qualitatively similar 
behavior. Thus, regardless of the strength of the four fold anisotropy field in an 
experimental sample, a change in the order of the switching transition should be 
observed as a function of temperature. 

Experimental evidence from an infrared spectroscopy (IR) study supports the 
scenario for a second order transition at room temperature. In their study, Arrecis 
et al. [25] observe only one peak in their IR spectra over the coverage range 
6 e [0.044,0.22]. If there were coexisting phases present, this would show as 
additional peaks in the spectra corresponding to the difi'erent symmetries of the 
phases. The LEED study of Griffiths et al. was performed at lower temperatures 
(T ~ 200 K) although the temperature was not held constant [24]. Their data 
shows indications of a coexistence region between 9 G [0.05,0.16] which could be 
interpreted in favor of a first order transition. The temperature range in Ref. 24 
is somewhat higher than our predictions but is in qualitative agreement with our 
work. Very recently, Okwamoto [43] has presented results of mean field calculations 
on the model of Lau and Ying [27] showing the switching transition to be of first 
order at all temperatures below the order - disorder line. However, it should be 
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noted that as the cubic anisotropy field vanishes at the transition, Okwamoto's 
method is equivalent to treating the pure XY model with a mean field theory. In 
this case it is well known that the mean field treatment does not predict the order 
of the transition correctly as the angular spin fluctuations destroy conventional 
long range order. 

The major drawbacks in this study were the slow dynamics of the model and the 
finite size effects. These together made it impossible for us within given computer 
time to locate the multicritical point more accurately. In addition, the adsorbate 
induced part generates all allowed anisotropy fields which complicates the analysis. 
Prom an experimental point of view, however, the location of the multicritical point 
may vary from sample to sample, depending on the possible intrinsic eight fold 
field. Thus a more accurate theoretical determination of may not prove to be 
necessary after all. 
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Appendix 

In this Appendix, we will discuss the hydrogen induced anisotropy fields in more 
detail. We shall first write a series expansion for the part of the Hamiltonian which 
is due to the adsorbate. We shall then expand the H — W interaction potential in 
terms of the anisotropy fields, and discuss the convergence of this expansion. 

We start from the adsorbate induced part of Eq. (7) and denote fugacity by 
7 = exp The summation i goes over all lattice sites and v over both bridge 
sites on each lattice site. Notation is in all cases equivalent to that in the text. 



jjeff _ jjclean ^ _/j^T^log(l + e''^^^'''^-''^) (A.l) 

i,v 
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Next, we study the effect of the H — W interaction terms which 
generate all the anisotropy fields induced by the hydrogen degrees of freedom. 
Due to symmetry arguments, the odd powers of u vanish. Let us now for the sake 
of argument consider the term in Eq {A.5) for n = 4: 
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Y.^y^,^)' = E [E («^^. ■ %)E K^/^ • 

ijV ijf 3 k 

E ("i^m • ^m)E (^'i^^n " ^n)] , (^-7) 

m n 

which reduces to 

when we consider only the on-site terms. 

On the other hand, based on symmetry arguments we can readily write the 
expansion for (^.1) in a general form when we consider contributions only from 
the on-site terms: 

CO 

geff _ gclean = ^ ^ /i«^(T, ^) COs(4A;(/)0 (A.9) 
i k=0 

where (pi is the displacement angle. The terms where k — 1 and 2 correspond to 
the analysis in Sec. 2.3. 

When we now expand Eq. {A.8) in terms of the displacement angle (p, and note 
that u = {uxiUyY' , \u\ = uo, Ux = ttocos0, Uy = uosmcp and \oii^j\ = a, we get a 
contribution for the four fold anisotropy field as 



However, it is important to notice that this is not the total adsorbate induced 
field h'l{T,ii). Namely, the evaluation of higher powers of {Vi^y)^ in Eq. (A. 5) 
results in contributions not only to /j.^, but also to lower order anisotropy fields. 
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These contributions can be of the same order of magnitude as the leading terms 
which severely hampers the analysis of these expansions. Thus, for example, for 
the coefficient of the four fold anisotropy field we obtain 



kBT(auo\^ 3kBTfauo\^ 



1211024 



where Fn{n,T) is defined by {A.5) and (A.Q). 

To study analytically the interplay between the cubic and eighth order anisotropy 
fields, and the vanishing of the total /i4(T, //) at the switching transition, one 
must first ensure the overall convergence condition exp{— /?(Vj^i^ — fi)} < 1. 
Unfortunately, for the values we have used for the parameters a,/! and T, this 
condition is not usually met in the vicinity of the switching transition line. We 
have also verified this numerically. In addition, for the parameters that we use, 
the term auo/{kBT) is never less than one for physically reasonable temperatures, 
thus further hampering the convergence of (A. 11). The additional effect of the 
more complicated off-site terms, which involve products of cosines at different 
lattice sites, is also difficult to determine accurately but it seems evident that they 
should be included in a quantitatively accurate calculation. 
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Figure captions 
Figure 1 

Here we show schematically the relative positions of W and H atoms on the 
Vr(lOO) surface. The solid dots (•) denote the ideal W atom positions, and the 
vectors ttj denote the positions of the displaced W atoms. The hydrogen atoms 
reside on the bridge sites (x) i' = (i, z^) where i denotes the tungsten atom lattice 
sites, and u = 1,2 the two bridge sites related to each lattice site. The H atoms 
feel the W atom potential of the nearest sites only, e.g. Vn = —auix + aukx, 
cf. Eq. (7). We have also depicted the magnitudes of the displacements and the 
lattice constant. 

Figure 2 

(a). The eight directions favored by the positive and negative four fold anisotropy 
fields. The dotted directions are favored by a negative /i4 and the solid lines show 
directions favorable for a positive h4. A negative eight fold field hg favors all 
these directions, (b). The eight directions favored by a positive eight fold field. 
Compared to (a) these directions are rotated by |-. 

Figure 3 

(a). MCRG cumulant behavior typical for a second order transition at T = 1.4, 
and (b) for a first order behavior at low temperatures (T = 0.5) for parameter set 
I. See text for details. 

Figure 4 

Coverage (Eq. (8)) increases as a function of the chemical potential. At elevated 
temperatures the increment takes place smoothly whereas at lower temperatures 
an abrupt increase is seen. The latter would indicate a possible first order 
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transition at low temperatures as discussed in the text. The calculation is 
performed for parameter set I. 

Figure 5 

The phase diagrams for H/W (100) in the (T, //) plane, based on numerical 
simulations of our model Hamiltonian for parameter sets I (a) and II (b) The 
solid dots denote the MCRG results. The switching transition changes order at 
Tf above which the transition is of second order (solid line) , and below which it is 
assumed to be of first order (dashed line). The dotted line shows the mean field 
result. The solid line between the ordered and disordered phases in (a) is only a 
guide to the eye. 
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